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ABSTRACT 

Studies of transiting extrasolar planets offer an unique opportunity to get to know the internal 
structure of those worlds. The transiting exoplanet XO-5 b was found to have an anomalously high 
Safronov number and surface gravity. Our aim was to refine parameters of this intriguing system 
and search for signs of transit timing variations. We gathered high-precision light curves of two 
transits of XO-5 b. Assuming three different limb darkening laws, we found the best-fitting model 
and redetermined parameters of the system, including planet-to-star radius ratio, impact parameter 
and central time of transits. Error estimates were derived by the prayer bead method and Monte 
Carlo simulations. Although system's parameters obtained by us were found to agree with previous 
studies within one sigma, the planet was found to be notable smaller with the radius of 1.03+qqj 
Jupiter radii. Our results confirm the high Safronov number and surface gravity of the planet. With 
two new mid-transit times, the ephemeris was refined to BJD TDB = (2454485. 66842± 0.00028) + 
(4. 1877537 ± 0.000017)£ . No significant transit timing variation was detected. 
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1. Introduction 

Transiting exoplanets are of great importance for the astrophysics of extrasolar 
planetary systems. If the inclination of the planetary orbit is close to 90° , a planet 
periodically moves across a stellar disk, blocking a small fraction of the flux. This 
phenomenon - a transit - is observed as a small (up to ~3%) drop of host star's 



'Based on observations collected at the Centra Astronomico Hispano Aleman (CAHA), oper- 
ated jointly by the Max-Planck Institut fiir Astronomie and the Instituto de Astrofisica de Andalucia 
(CSIC). 
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brightness. Combining spectroscopic and photometric data gives unique oppor- 
tunity to determine planetary mass and radius, and hence the surface gravity and 
average density - a key parameter for studying the internal structure of exoplanets. 
The first exoplanetary transits were observed for HD 209458 b (Charbonneau et al. 
2000; Henry at al. 2000) which was discovered with the radial-velocity technique 
(Mazeh et al. 2000; Henry at al. 2000). The efficiency of photometric searches 
for transiting planets is seriously affected by the high fraction of false positive sce- 
narios because variety of phenomena may mimic a planetary transit light curve. 
Examples of these false positives are central transit of a low-mass star in front of a 
large main-sequence star or red giant, grazing eclipses in systems comprising two 
main-sequence stars or a contamination of a fainter eclipsing binary along the same 
line of sight (Charbonneau et al. 2004). OGLE-TR-56 b is the first exoplanet dis- 
covered by the transit method. It was initially listed as a planetary candidate found 
by the Optical Gravitational Lensing Experiment (OGLE, Udalski et al. 1997) in 
the direction of the Galactic Centre (Udalski et al. 2002a,b). Then, the plane- 
tary nature of OGLE-TR-56 b was confirmed by radial velocity measurements by 
Konacki et al. (2003). Since then, over 100 extrasolar transiting planets have been 
discovered by numerous surveys. 

High-precision photometric follow-ups of transiting exoplanets allow to refine 
planetary and stellar parameters. Such observations may also lead to discoveries of 
additional, even very low mass bodies in extrasolar systems if transit time variations 
(TTVs) are detected (Miralda-Escude 2002; Holman & Murray 2005; Agol et al. 
2005). The example of the Kepler 9 planetary system, in which there are at least 2 
transiting planets, clearly demonstrates the usefulness of the TTV method (Holman 
et al. 2010). Preliminary detections of the TTV signal were reported for the WASP- 
3 b and WASP-10 b transiting planets (Maciejewski at al. 2010a,b). The source 
of these deviations from a strictly Keplerian case could be additional planets in 
these planetary systems. Detection of a TTV signal together with transit duration 
variations (TDVs) shifted in a phase by 7l/2 would indicate the presence of an 
exomoon of the transiting planet (Kipping 2009). 

The star XO-5 (GSC 02959-00729, a = 07 h 46 m 52 s , 5 = ±39°05'41") was 
found to harbour a transiting hot Jupiter, XO-5 b, which orbits its host star within 
4.2 days and whose mass and radius were found to be Mb = 1.15 ±0.08 Mj and 
Rb = 1.15 ±0.12 Rj, respectively (Burke et al. 2008). High-resolution spectral 
observations revealed that the host star is a dwarf of the G8 spectral type, has the 
effective temperature of 5510 ±44 K and is located 270 ±25 pc from the Sun. 
The age of the system was estimated to be 8.5 ±0.8 Gyr. The planet was inde- 
pendently confirmed by Pal et al. (2009) who refined system's parameters. The 
planetary mass and radius were found to be 1.059 ±0.028 Mj and 1.109 ±0.050 
Rj , respectively. The host star was found to be slightly cooler with the effective 
temperature of 5370 ±70 K and more evolved with the age of 14.8 ±2.0 Gyr. The 
orbital period was refined to 4.187757 ±0.000011 d and no sign of variations in 
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transit timing was found. Using theoretical models, Pal et al. (2009) found that the 
planet's core could be smaller than 10 Earth masses if a planet older than 4.5 Gyr 
is considered. 

In this paper we present results of high-precision photometric follow-up obser- 
vations whose aim was to refine parameters of the XO-5 system. 

2. Observations and data reduction 

Two transits of XO-5 b were observed with the 2.2-m telescope at Calar Alto 
Observatory (Spain) during 2 runs on 2010 January 21 and November 10. The 
Calar Alto Faint Object Spectrograph (CAFOS) in imaging mode was used as a 
detector. It was equipped with the SITe CCD matrix (2048 x 2048, 24,um pixel, 
0.53 arcsec per pixel). A subframe limiting the field of view to 5.3 x 3.7 arcmin 
was used to shorten the read-out time. During the second run, binning in 2 x 2 
mode was used to additionally shorten the read-out time. The subframe was chosen 
in such a way to observe simultaneously a nearby comparison star GSC 02959- 
01873 which has a colour index and brightness similar to the ones of XO-5. This 
criterion was expected to minimise photometric trends caused by the differential 
atmospheric extinction. The photometric monitoring was performed in the John- 
son /?-band filter in which the instrument set-up is most sensitive. The telescope 
was significantly defocused and hence stellar profiles exhibited a donut-like shape 
(Fig. 1). This method minimises random and flat-fielding errors (e.g. Southworth 
et al. 2009). The stellar flux was spread over a ring of ~ 17 arcsec in diameter. 
A visual inspection of the Digitized Sky Survey (DSS) images revealed no faint 
neighbour stars up to 40 arcsec in the direct vicinity of XO-5 nor the comparison 
star down to the limiting magnitude which is at least R ~ 20.5 mag. The stellar 
images were kept exactly at the same position in the CCD matrix during each run 
thanks to auto guiding. Precise timing was assured by synchronising the computer's 
clock to Coordinated Universal Time (UTC) by Network Time Protocol software, 
accurate to better than 0.1 s. 

CCD frames were processed using a standard procedure including debiassing 
but not flat-fielding which does not improve measurements in the case of defo- 
cused and auto-guided observations but may also degrade the quality of data (see 
e.g. Southworth et al. 2010). Indeed, our tests showed that dividing by a flat-field 
frame had no detectable effect on a final light curve. The magnitudes of XO-5 
and comparison star were determined with differential aperture photometry. The 
aperture radii ranged from 8 to 18 pixels and 6 to 14 pixels for runs 1 and 2, re- 
spectively. Beginning from the aperture radii of 8 and 12 pixels for runs 1 and 2, 
respectively, the photometric scatter exhibited a flat minimum (or rather a plateau), 
in practice insensitive to greater aperture sizes. To avoid possible neighbour-star 
contamination, the light curve with minimal apertures for which the lowest scatter 
was achieved (i.e. 8 and 12 pixels for runs 1 and 2, respectively), were taken as the 
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Fig. 1. The flux distribution of the stellar profile. The donut-like shape is caused by defocusing the 
telescope. The 30 x 30 pixel fragment corresponds to 15.9 x 15.9 arcsec on the sky (no binning). 
The logarithmic vertical scale was used for better contrast. The data were extracted from a typical 
50-s exposure collected during the run on 2010 January 21. 

final ones. The comparison star was used to set the zero level of the photometric 
scale. No trends which could be approximated by a first- or second-order polyno- 
mial were detected in both transit light curves due to small changes of the airmass. 
Therefore, no detrending procedure was applied. 

During the first run, sky conditions were non-photometric and data were ac- 
quired in gaps between clouds. The ingress phase was lost and egress data were 
significantly affected by thin clouds. Therefore, these data were found to be of 
marginal use. In the second run, almost photometric conditions allowed to record a 
complete transit light curve. Occasionally thin clouds affected flat-bottom phase of 
the transit. The exposure time was set to get roughly 1-min cadence and avoid satu- 
ration. It was refined during out-of-transit phase to achieve maximal efficiency. To 
avoid effecting transit timing, no exposure-time changes were done during ingress 
nor egress phases. Observations in both runs were carried out in dark time, just a 
few days before or after a new moon. Hence they were not affected by moonlight. 
The details of both observing runs are presented in Table 1 . 



3. Results 



Individual light curves were modelled with the JKTEBOP code (Southworth et 
al. 2004a, 2004b) which is based on the EBOP programme (Eclipsing Binary Orbit 
Program; Etzel 1981; Popper & Etzel 1981). The software models both compo- 
nents of a system - a planet and a host star - as biaxial ellipsoids and performs a 
numerical integration in concentric annuli over the surface of each body to obtain 
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Table 1 

The summary of observing runs: 7V exp - the number of useful exposures, X - airmass changes 
during a given run, r exp - exposure times. Dates are given in UT at the beginning of nights. 



Run 


Date 


A'exp 


X 




?exp (s) 


Binning mode 


1 


2010 January 21 


143 


1.07-> 1.00 - 


+ 1.18 


45, 50, 60 


lxl 


2 


2010 November 10 


190 


1.16 ->• 1.00- 


> 1.02 


40, 45, 50 


2x2 



the flux coming from the system. This feature of the code allows to avoid small and 
spherical planet approximations which are used in analytic light-curve generators 
based on Mandel & Agol (2002), and hence to derive planet's oblateness. A model 
is fitted to the data by the Levenberg-Marquardt least-square procedure. The code 
converges rapidly toward a reliable solution and diminishes the correlation between 
fitted parameters (Southworth 2008). 

For the high-quality light curve from the second run, five parameters describ- 
ing a shape of a light curve were allowed to float during fitting procedure. We 
used fractional radii of the host star and planet, defined as r» = y and h> = ^ , 
respectively, where /?* and R\, are the absolute radii of the bodies and a is the or- 
bital semi-major axis. In practise, the combinations of these parameters were used: 
a sum r* + and ratio k = r*/r\, because they were found to be the less corre- 
lated with each other (Southworth 2008). The directly fitted orbital inclination, i, 
allowed to calculate the transit parameter b = cos i . The initial values of parame- 
ters listed above were taken from Pal et al. (2009). To detect any variation in transit 
timing between observed transits, the mid-transit time was set as a free parameter 
whose initial value was calculated according to the ephemeris given by Burke et al. 
(2008). High-quality light curves need to set limb-darkening coefficients (LDCs) 
as free parameters and test different limb-darkening (LD) laws (Southworth 2008). 
Hence, we considered linear, logarithmic and square-root LD, for which theoreti- 
cal LDCs in the Johnson R band were bilinearly interpolated from tables by Van 
Hamme (1993). In a first iteration we allowed both linear u and non-linear v co- 
efficients to vary but unrealistic values were derived. Therefore, we kept v fixed at 
theoretical values and allow only u to vary in final fitting runs. The contribution of 
v into the error budget was included by perturbing it by ±0.1 around its theoretical 
value and assuming a flat distribution. 

The errors of derived parameters were determined in two ways for each combi- 
nation of the data set and adopted LD law. Firstly, we run 1000 Monte Carlo (MC) 
simulations and a spread range of a given parameter within 68.3% was taken as its 
error estimate. Secondly, the prayer-bead method (e.g. Desert et al. 2009; Winn 
et al. 2009) was used to check whether red noise is present in our data. MC errors 
were found to be 2-3 times smaller than the values returned by the prayer bead 
method, thus the latter ones were used as the final values. This finding indicates 
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that beside Poisson noise there is correlated (red) noise in the light curve. 
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Fig. 2. Light curves of two transits of XO-5 b observed on 2010 January 21 and November 10. The 
best-fitting models, based on the square-root limb darkening law, are plotted with continuous lines. 
The residuals are shown in bottom panels. 



Table2 

Parameters of transit light-curve modelling derived for the high-quality light curve acquired in the 
second run on 2010 November 10. Linear, logarithmic and square-root limb darkening laws were 
considered. Tq is based on UTC and is given as 245551 1 + . See text for details. 



Parameter 


Linear 


Logarithmic 


Square-root 


r* + r b 


Uil 11 -0.0036 




1102+ UUU3y 


k 


n ims+ 00014 
u.iui»_ 00014 


0-1018jffi 


1 013+ 00014 


r* 


1009+ 0035 


1 032+ 0037 


1 000+ 00034 
u.iuuu_ 00040 


fb 


n nin?s+ 00052 
0.01028_ 000045 


01050+ 00053 


0.01013^;K 


i (deg) 


86.801°;* 


86.61^ 


86.9+0;* 


b{R*) 


53+ 08 
u - :):, -0.08 


5^+ 09 
u - 3:, -0.09 


51+ 08 


u 


45+ 07 


57+ 08 
"•-"-0.09 


12+ 007 


V 




0.2171 a 


0.5419 a 


Tq (JDijtc) 


0.66486_ 00030 


0.66487_ 00033 


0.66483_ 000029 


rms (mmag) 


0.8334 


0.8336 


0.8329 


Xred 


1.6325 


1.6326 


1.6304 



a permuted by ±0. 1 on a flat distribution 



Table 2 contains results obtained for the second-run light curve and individual 
LD laws. The sub-millimagnitude precision was achieved with the rms of 0.83 
mmag. The smallest %^ ed was achieved for the square-root LD law and this set of 
system's parameters was taken as the final one. However, it is worth noting that 
values of rms of individual fits are comparable to each other what suggests that the 
choice of the LD law is not significant for the analysed light curve. Moreover, the 
timing errors seem to be independent on the choice of the LD law. In all three cases 
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the precision between 20 and 29 s was achieved. 

Full analysis of the first-run light curve was found to be unreliable, thus all 
parameters but mid-transit time were fixed and this light curve was used only for 
transit timing. The values of system's parameters were adopted from best-fitting 
model obtained for the second-run light curve. The error of mid-transit time was 
determined with the prayer bead method. The %j 2 ed of the fit was found to be 3.8 
with the rms =1.4 mmag what reflects the lower quality of the data. 

The light curves acquired in both runs are plotted in Fig. 2 together with the 
best-fitting models and residual^. 

The mid-transit times available in the literature and 2 new determinations re- 
ported in this paper were used to refine the transit ephemeris. We used 8 high- 
quality mid-transit times from Burke et al. (2008) rejecting unreliable ones which 
are based on the original XO photomeuy. We also used 4 mid-transit times given 
by Pal et al. (2009); values for the individually fitted transit centres were taken. 
The mid-transit times were transformed from JD based on UTC into BJD based 
on Barycentric Dynamical Time (TDB) using the on-line converted by Eastman et 
al. (2010). As a result of fitting a linear function of the epoch and period P^, we 
obtained: 

T = 2454485.66842 ± 0.00028 BJD TD b 

P h = 4. 1877537 ±0.0000017 d. 

Individual mid-transit errors were taken as weights. Results for new mid-transit 
times are summarised in Table 3. The observation minus calculation (O — C) di- 
agram, plotted in Fig. 3, shows no significant deviation of data points from the 
new linear ephemeris. Some literature measurements deviate by more than 1 sigma 
which may be a result of underestimated timing errors. 



Table3 

Results of transit timing. To denotes the mid-transit times given as JD (based on Coordinated 

Universal Time, UTC) and BJD (based on Barycentric Dynamical Time, TDB). Errors of 
mid-transit times are in days. O — C values were calculated according to the new ephemeris. 



Run 


T (JDutc) 


T (BJD TDB ) 


To error 


Epoch 


O-C(d) 


1 

2 


2455218.5194 
2455511.66483 


2455218.5255 
2455511.66809 


+0.0013 
-0.0008 
+0.00023 
-0.00029 


175 
245 


+0.0002 
+0.00001 



2 The data in a machine-readable form are available at http://web.astri.umk.pl/ttv 
3 http://astroutils. astronomy.ohio-state.edu/time/utc 2bjd.html 
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Fig. 3. The observation minus calculation (O — C) diagram for transit timing of XO-5 b, generated 
according to the new linear ephemeris. The open triangles and squares denote individual mid-transit 
times published by Burke et al. (2008) and Pal et al. (2009), respectively. The filled symbols mark 
two mid-transit times reported in this paper. No transit timing variation was detected. 

4. Physical properties of XO-5 system 

Results of light-curve modelling allowed us to calculate planetary, stellar and 
geometrical parameters. We used results from the best-fitting model based on 
the second-run data and square-root LD law. R\, was calculated assuming a = 
0.0488 ±0.0006 AU (Pal et al. 2009). The planetary mean density p b uses XO- 
5 b's mass of 1.059 ±0.028 Mj, known from radial velocity measurements (Pal 
et al. 2009). Surface gravitational acceleration, g\, , was calculated according to a 
formula (Southworth et al. 2007): 

2n Vl-e\ 
8b = tt o . . K*, (1) 
Pb rgsmz 

where F\> is the orbital period of the planet (see Sect. 3) and K* is the stellar velocity 
amplitude taken from Pal et al. (2009). The orbital eccentricity was assumed to be 
zero (Burke et al. 2008; Pal et al. 2009). 

The equilibrium temperature, T sq , may be derived assuming the effective tem- 
perature of the host star derived by Pal et al. (2009) and using the relation (South- 
worth 2010): 

where A and F are the Bond albedo and heat redistribution factor, respectively. 
This formula was simplified by assuming relation A = 1 — 4F and the modified 
equilibrium temperature r e ' q was calculated (Southworth 2010). 

We calculated the Safronov number which determines the efficiency with 
which a planet gravitationally scatters other bodies (Safronov 1972). This parame- 
ter is proportional to the ratio of the escape velocity from the planet and velocity of 
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© = ttV , (3) 



its orbital motion and can be calculated according to a formula (Southworth 2010): 

where the the value of 0.88 ± 0.03 M was taken as star's mass M* (Pal et al. 
2009). 

Assuming that the planet's rotation and orbital periods are synchronised, the 
oblateness of the planet, defined as: 

/=!--, (4) 

n 

where r\ and r2 are the polar and equatorial radii, respectively, was found to be 
0.0014. Such a small value indicates that any deviations of the XO-5 b's shape 
from an ideal sphere may be neglected. 

Results of calculations are collected in Table 4 in which literature determina- 
tions are also given for comparison. 



Table4 

Physical properties of the XO-5 system derived from light-curve modelling. Listed parameters are 
explained in the text. Values derived by Burke et al. (2008) and Pal et al. (2009) are given for 

comparison. 



Parameter 


This work 


Burke et al. (2008) 


Pal et al. (2009) 




Planetary properties 




*b (Ri) 


1 03+ 006 


1.15±0.12 


1.109 ±0.050 


p b (g cm" 3 ) 
gb (m s~ 2 ) 


1 97+O.19 
1 - z/ -0.17 

24.6+Jj 


1.02±0.3 
22±5 


96+°- 14 
ayD -0.11 

21 4+ 2 ' 1 


Fe'q (K) 


1201+^ 


1244 ±48 


1221 ±27 




/ 


114+ 0013 
U- 1 ^-0.01 1 

0.0014 


0.10±0.01 


0.105 ±0.005 




Geometrical parameters 




( (deg) 


86.9 ±0.4 


86.8±0.9 


86.7 ±0.4 


b(R*) 
Rb/R* 


0.51 ±0.08 
1013+ 00014 


0.55 ±0.09 
0.106 ±0.003 


5 62+ 033 
0.1050 ±0.0009 




Stellar properties 




R* {R&) 


1 05+ 005 
1 - u; '-0.04 


1.11 ±0.09 


1.08 ±0.04 


P* (Po) 


76+ 07 


0.72 ±0.14 




log#» (cgs) 


4 34+O.O6 
^■-^-0.05 


4.34 ±0.07" 


4.31 ±0.03 



a Here spectroscopic determination was taken. 



5. Concluding discussion 

High-precision photometric follow-up observations of transiting planets are 
crucial for verifying physical properties of their planetary systems. Although sys- 
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tern's parameters obtained by us were found to agree with previous studies by Burke 
et al. (2008) and Pal et al. (2009) within one sigma, the planet was found to be no- 
table smaller, i.e. by 7%, comparing to the value obtained by Pal et al. (2009). A 
shallower transit, and hence smaller ratio of planet and star radii, may be a sign of 
the presence of a third light contributing to the system's light curve, e.g. a faint 
star blended with XO-5 in the DSS image. To check this scenario, we explored the 
Naval Observatory Merged Astrometric Dataset (NOMAD, Zacharias et al. 2004) 
and found a faint source (R = 19.62 mag), located 9 arcsec to the south from XO-5. 
Assuming an extreme and unlikely case, in which a total flux of the neighbour ob- 
ject contributed to the XO-5 light curve, the planetary radius would be affected by 
only 0.05% - over 100 times less than the detected difference. This result strength- 
ens the reliability of our finding. 

The smaller planet's radius induces larger values of the mean density, pb = 
1.27^Q'}y g cm~ 3 , and gravitational acceleration, g\, = 24.6^' 3 m s 2 > whose 
values are comparable to Jupiter's ones. Our results confirm the high Safronov 
number of XO-5 b which is equal to 0.114+qqJj . Such a high value suggests that 
the planet has not undergone a phase of the significant mass evaporation driven by 
the radiation of the host star (e.g. Hansen & Barman 2007). This finding is not 
surprising if one considers the spectral type of the star (G8) and a relatively small 
amount of the extreme ultraviolet radiation that it emits. 

The transit times were found to be periodic within a 3-sigma level. This finding 
reveals no hint for the existence of the additional planet which could be discovered 
by the TTV method in the XO-5 system. 
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